Detection of Cytoplasmic and Nuclear Circular RNA via RT-qPCR.

Circular RNA (circRNA) is an intriguing class of non-coding RNA that exists as a continuous closed loop. With the improvements in high throughput sequencing, biochemical analysis, and bioinformatic algorithms, studies on circRNA expression became abundant in recent years. However, functional studies of circRNA are still limited. Subcellular localization of circRNA may provide some clues in elucidating its biological functions by performing subcellular fractionation assay. Notably, circRNAs that are predominantly found in the cytoplasm are more likely to be involved in post-transcriptional gene regulation, e.g., acting as micoRNA sponge, whereas nuclear-retained circRNAs are predicted to play a role in transcriptional regulation. Subcellular fractionation could help researchers to narrow down and prioritize downstream experiments. The majority of the currently available protocols describe the steps for subcellular fractionation followed by western blot analysis for protein molecules. Here, we present a protocol for the subcellular fractionation of cells to detect circRNA via RT-qPCR with divergent primers. Moreover, detailed steps for the generation of specific circRNAs-enriched cDNA included in this protocol will enhance the amplification and detection of low-abundance circRNAs. This will be useful for researchers studying low-abundance circRNAs. Key features This protocol builds upon the method developed by Gagnon et al. (2014) and extends its application to circRNA study. Protocol for amplification of low levels of circRNA expression. Analysis takes into consideration the ratio of cytoplasmic RNA concentration to nuclear RNA concentration.

Bat ASC2 suppresses inflammasomes and ameliorates inflammatory diseases.

Bats are special in their ability to live long and host many emerging viruses. Our previous studies showed that bats have altered inflammasomes, which are central players in aging and infection. However, the role of inflammasome signaling in combating inflammatory diseases remains poorly understood. Here, we report bat ASC2 as a potent negative regulator of inflammasomes. Bat ASC2 is highly expressed at both the mRNA and protein levels and is highly potent in inhibiting human and mouse inflammasomes. Transgenic expression of bat ASC2 in mice reduced the severity of peritonitis induced by gout crystals and ASC particles. Bat ASC2 also dampened inflammation induced by multiple viruses and reduced mortality of influenza A virus infection. Importantly, it also suppressed SARS-CoV-2-immune-complex-induced inflammasome activation. Four key residues were identified for the gain of function of bat ASC2. Our results demonstrate that bat ASC2 is an important negative regulator of inflammasomes with therapeutic potential in inflammatory diseases.

Targeting the inflammasome in Parkinson's disease.

Parkinson's disease (PD) is one of the most common neurodegenerative diseases in which neuroinflammation plays pivotal roles. An important mechanism of neuroinflammation is the NLRP3 inflammasome activation that has been implicated in PD pathogenesis. In this perspective, we will discuss the relationship of some key PD-associated proteins including α-synuclein and Parkin and their contribution to inflammasome activation. We will also review promising inhibitors of NLRP3 inflammasome pathway that have potential as novel PD therapeutics. Finally, we will provide a summary of current and potential in vitro and in vivo models that are available for therapeutic discovery and development.

Out of chaos: Phylogenomics of Asian Sonerileae.

Sonerileae is a diverse Melastomataceae lineage comprising ca. 1000 species in 44 genera, with >70% of genera and species distributed in Asia. Asian Sonerileae are taxonomically intractable with obscure generic circumscriptions. The backbone phylogeny of this group remains poorly resolved, possibly due to complexity caused by rapid species radiation in early and middle Miocene, which hampers further systematic study. Here, we used genome resequencing data to reconstruct the phylogeny of Asian Sonerileae. Three parallel datasets, viz. single-copy ortholog (SCO), genomic SNPs, and whole plastome, were assembled from genome resequencing data of 205 species for this purpose. Based on these genome-scale data, we provided the first well resolved phylogeny of Asian Sonerileae, with 34 major clades identified and 74% of the interclade relationships consistently resolved by both SCO and genomic data. Meanwhile, widespread phylogenetic discordance was detected among SCO gene trees as well as species trees reconstructed using different tree estimation methods (concatenation/site-based coalescent method/summary method) or different datasets (SCO/genomic/plastome). We explored sources of discordance using multiple approaches and found that the observed discordance in Asian Sonerileae was mainly caused by a combination of biased distribution of missing data, random noise from uninformative genes, incomplete lineage sorting, and hybridization/introgression. Exploration of these sources can enable us to generate hypotheses for future testing, which is the first step towards understanding the evolution of Asian Sonerileae. We also detected high levels of homoplasy for some characters traditionally used in taxonomy, which explains current chaotic generic delimitations. The backbone phylogeny of Asian Sonerileae revealed in this study offers a solid basis for future taxonomic revision at the generic level.

Expression Level Dominance and Homeolog Expression Bias Upon Cold Stress in the F1 Hybrid Between the Invasive Sphagneticola trilobata and the Native S. calendulacea in South China, and Implications for Its Invasiveness.

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid's responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

Glacial vicariance and oceanic circulation shape population structure of the coastal legume Derris trifoliata in the Indo-West Pacific.

PREMISE: The phylogeography of coastal plant species is shaped by contemporary and historical biogeographic processes. In this study, we aim to decipher the phylogeography of Derris trifoliata, a woody legume of relatively recent origin and wide distribution, in coastal areas in the Indo-West Pacific (IWP) region. METHODS: Genetic diversity and population structure were assessed by analyzing six nuclear and three chloroplast DNA sequences from 30 populations across the species' range. Phylogeography was inferred by estimating gene flow, divergence time, historical population size changes, and historical habitat suitability using paleoclimatic niche modeling. RESULTS: High genetic diversity was observed at the species level. The populations of three oceanic regions included in this study (i.e., Indian Ocean, South China Sea, and Pacific Ocean) formed distinct clades and likely diverged during the late Pleistocene. Potential barriers to gene flow were identified, including the Sunda and Sahul shelves, geographic distance, and current patterns of oceanic circulation. Analysis of changes in population size supported the bottleneck model, which was strengthened by estimates of habitat suitability across paleoclimatic conditions. CONCLUSIONS: The once widespread distribution of D. trifoliata was fragmented by changes in climatic suitability and biogeographic barriers that arose following sea-level changes during the Pleistocene. In addition, contemporary patterns of oceanic circulation and geographic distance between populations appear to maintain genetic differentiation across its distribution in the IWP.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea.

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

Contrasting Phylogeographic Patterns in Lumnitzera Mangroves Across the Indo-West Pacific.

Mangroves are ecologically important forest communities in tropical and subtropical coasts, the effective management of which requires understanding of their phylogeographic patterns. However, these patterns often vary among different species, even among ecologically similar taxa or congeneric species. Here, we investigated the levels and patterns of genetic variation within Lumnitzera consisting of two species (L. racemosa and L. littorea) with nearly sympatric ranges across the Indo-West Pacific (IWP) region by sequencing three chloroplast DNA regions (for both species) and genotyping 11 nuclear microsatellite loci (for L. littorea). Consistent with findings in studies on other mangrove species, we found that both L. racemosa and L. littorea showed relatively high genetic variation among populations but low genetic variation within populations. Haplotype network and genetic clustering analyses indicated two well-differentiated clades in both L. racemosa and L. littorea. The relationship between geographic and genetic distances and divergence time estimates of the haplotypes indicated that limited dispersal ability of the propagules, emergence of land barriers during ancient sea-level changes, and contemporary oceanic circulation pattern in the IWP influenced the current population structure of the two species. However, the position of genetic break was found to vary between the two species: in L. racemosa, strong divergence was observed between populations from the Indian Ocean and the Pacific Ocean possibly due to land barrier effect of the Malay Peninsula; in L. littorea, the phylogeographic pattern was created by a more eastward genetic break along the biogeographic barrier identified as the Huxley's line. Overall, our findings strongly supported previous hypothesis of mangrove species divergence and revealed that the two Lumnitzera species have different phylogeographic patterns despite their close genetic relationship and similar current geographic distribution. The findings also provided references for the management of Lumnitzera mangroves, especially for the threatened L. littorea.

Identification and characterization of a novel Epstein-Barr Virus-encoded circular RNA from LMP-2 Gene.

Epstein-Barr virus (EBV) has been recently found to generate novel circular RNAs (circRNAs) through backsplicing. However, comprehensive catalogs of EBV circRNAs in other cell lines and their functional characterization are still lacking. In this study, we have identified a list of putative EBV circRNAs in GM12878, an EBV-transformed lymphoblastoid cell line, with a significant majority encoded from the EBV latent genes. A novel EBV circRNA derived from the exon 5 of LMP-2 gene which exhibited highest prevalence, was further validated using RNase R assay and Sanger sequencing. This circRNA, which we term circLMP-2_e5, can be universally detected in a panel of EBV-positive cell lines modelling different latency programs. It ranges from lower expression in nasopharyngeal carcinoma (NPC) cells to higher expression in B cells, and is localized to both the cytoplasm and the nucleus. We provide evidence that circLMP-2_e5 is expressed concomitantly with its cognate linear LMP-2 RNA upon EBV lytic reactivation, and may be produced as a result of exon skipping, with its circularization possibly occurring without the involvement of cis elements in the short flanking introns. Furthermore, we show that circLMP-2_e5 is not involved in regulating cell proliferation, host innate immune response, its linear parental transcripts, or EBV lytic reactivation. Taken together, our study expands the current repertoire of putative EBV circRNAs, broadens our understanding of the biology of EBV circRNAs, and lays the foundation for further investigation of their function in the EBV life cycle and disease development.

Comparison of eight complete chloroplast genomes of the endangered Aquilaria tree species (Thymelaeaceae) and their phylogenetic relationships.

Aquilaria tree species are naturally distributed in the Indomalesian region and are protected against over-exploitation. They produce a fragrant non-timber product of high economic value, agarwood. Ambiguous species delimitation and limited genetic information within Aquilaria are among the impediments to conservation efforts. In this study, we conducted comparative analysis on eight Aquilaria species complete chloroplast (cp) genomes, of which seven were newly sequenced using Illumina HiSeq X Ten platform followed by de novo assembly. Aquilaria cp genomes possess a typical quadripartite structure including gene order and genomic structure. The length of each of the cp genome is about 174 kbp and encoded between 89 and 92 proteins, 38 tRNAs, and 8 rRNAs, with 27 duplicated in the IR (inverted repeat) region. Besides, 832 repeats (forward, reverse, palindrome and complement repeats) and nine highly variable regions were also identified. The phylogenetic analysis suggests that the topology structure of Aquilaria cp genomes were well presented with strong support values based on the cp genomes data set and matches their geographic distribution pattern. In summary, the complete cp genomes will facilitate development of species-specific molecular tools to discriminate Aquilaria species and resolve the evolutionary relationships of members of the Thymelaeaceae family.

Molecular sexing and preliminary assessment of population sex ratio of the endangered Malayan tapir (Tapirus indicus) in Peninsular Malaysia.

A molecular sexing method by polymerase chain reaction (PCR) amplification of a portion of the sex-determining region Y (SRY) and the zinc finger (ZF) gene, as well as six equine Y-chromosome-specific microsatellite markers, were tested in the Malayan tapir (Tapirus indicus). While the microsatellite markers did not yield any male-specific amplicons for sex-typing, the SRY/ZF marker system produced reliable molecular sexing results by accurately sex-typing 31 reference Malayan tapirs, using whole blood, dried blood spot (DBS), or tissue samples as materials for DNA extraction. The marker system was also tested on 16 faecal samples, and the results were in general consistent with the pre-determined sexes of the animals, despite some amplification failures. A preliminary estimation of wild Malayan tapir population sex ratio was estimated from the Wildlife Genomic Resource Bank (WGRB) database of the Malaysian Department of Wildlife and National Parks (PERHILITAN), zoos, and the Sungai Dusun Wildlife Conservation Centre (WCC), as well as from the results of molecular sexing 12 samples of unknown sex. The overall sex ratio favoured females, but the deviation from parity was statistically not significant when tested using the binomial test (p > 0.05), which may be due to reduced statistical power caused by small sample sizes.

Cross-amplification of microsatellite markers across agarwood-producing species of the Aquilarieae tribe (Thymelaeaceae).

Tree species in the Aquilarieae tribe of the Thymelaeaceae family produce agarwood, a natural product highly valued for its fragrance, but the species are under threat due to indiscriminate harvesting. For conservation of these species, molecular techniques such as DNA profiling have been used. In this study, we assessed cross-amplification of microsatellite markers, initially developed for three Aquilaria species (A. crassna, A. malaccensis, and A. sinensis), on ten other agarwood-producing species, including members of Aquilaria (A. beccariana, A. hirta, A. microcarpa, A. rostrata, A. rugosa, A. subintegra, and A. yunnanensis) and Gyrinops (G. caudata, G. versteegii, and G. walla), both from the Aquilarieae tribe. Primers for 18 out of the 30 microsatellite markers successfully amplified bands of expected sizes in 1 sample each of at least 10 species. These were further used to genotype 74 individuals representing all the 13 studied species, yielding 13 cross-amplifiable markers, of which only 1 being polymorphic across all species. At each locus, the number of alleles ranged from 7 to 23, indicating a rather high variability. Four markers had relatively high species discrimination power. Our results demonstrated that genetic fingerprinting can be an effective tool in helping to manage agarwood genetic resources by potentially supporting the chain-of-custody of agarwood and its products in the market.

Analyses of Plastome Sequences Improve Phylogenetic Resolution and Provide New Insight Into the Evolutionary History of Asian Sonerileae/Dissochaeteae.

Sonerileae/Dissochaeteae (Melastomataceae) comprises ca. 50 genera, two thirds of which occur in Southeast Asia. Phylogenetic relationships within this clade remain largely unclear, which hampers our understanding of its origin, evolution, and biogeography. Here, we explored the use of chloroplast genomes in phylogenetic reconstruction of Sonerileae/Dissochaeteae, by sampling 138 species and 23 genera in this clade. A total of 151 complete plastid genomes were assembled for this study. Plastid genomic data provided better support for the backbone of the Sonerileae/Dissochaeteae phylogeny, and also for relationships among most closely related species, but failed to resolve the short internodes likely resulted from rapid radiation. Trees inferred from plastid genome and nrITS sequences were largely congruent regarding the major lineages of Sonerileae/Dissochaeteae. The present analyses recovered 15 major lineages well recognized in both nrITS and plastid phylogeny. Molecular dating and biogeographical analyses indicated a South American origin for Sonerileae/Dissochaeteae during late Eocene (stem age: 34.78 Mya). Two dispersal events from South America to the Old World were detected in late Eocene (33.96 Mya) and Mid Oligocene (28.33 Mya) respectively. The core Asian clade began to diversify around early Miocene in Indo-Burma and dispersed subsequently to Malesia and Sino-Japanese regions, possibly promoted by global temperature changes and East Asian monsoon activity. Our analyses supported previous hypothesis that Medinilla reached Madagascar by transoceanic dispersal in Miocene. In addition, generic limits of some genera concerned were discussed.

Comparative transcriptomics sheds light on differential adaptation and species diversification between two Melastoma species and their F1 hybrid.

Variation in gene expression has been shown to promote adaptive divergence, and can lead to speciation. The plant genus Melastoma, thought to have diversified through adaptive radiation, provides an excellent model for the study of gene expressional changes during adaptive differentiation and following interspecific hybridization. In this study, we performed RNA-seq on M. candidum, M. sanguineum and their F1 hybrid, to investigate the role of gene expression in species diversification within the genus. Reference transcriptomes were assembled using combined data from both parental species, resulting in 50 519 and 48 120 transcripts for the leaf and flower petal, after removing redundancy. Differential expression analysis uncovered 3793 and 2116 differentially expressed (DE) transcripts, most of which are between M. candidum and M. sanguineum. Differential expression was observed for genes related to light responses, as well as genes that regulate the development of leaf trichomes, a trait that among others is thought to protect plants against sunlight, suggesting the differential adaptation of the species to sunlight intensity. The analysis of positively selected genes between the two species also revealed possible differential adaptation to other abiotic stresses such as drought and temperature. In the hybrid, almost all possible modes of expression were observed at the DE transcripts, although at most transcripts, the expression levels were similar to that of either parent instead of being intermediate. A small number of transgressively expressed transcripts that matched genes known to promote plant growth and adaptation to stresses in new environments were also found, possibly explaining the vigour observed in the hybrid. The findings in this study provided insights into the role of gene expression in the diversification of Melastoma, which we believe is an important example for more cross-taxa comparisons in the future.

The complete chloroplast genome of the lipstick tree, Bixa Orellana (Bixaceae).

Bixa orellana is a small tree known for its red, oil-soluble pigment contained in the seed coat that is used as a natural dye and food coloring. In this study, we assembled and characterized the complete chloroplast genome of B. orellana as a resource for future genetic studies. With a total length of 159,825 bp, the chloroplast genome comprised of a large single-copy (LSC) region of 89,476 bp, a small single-copy (SSC) region of 19,617 bp, and two inverted repeat (IR) regions of 25,356 bp each. A total of 127 genes were predicted, consisting of 83 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. Phylogenetic analysis confirmed the position of B. orellana within the order Malvales.

The complete chloroplast genome sequence of Chengal (Neobalanocarpus heimii, Dipterocarpaceae), a durable tropical hardwood.

Known for its durable timber quality, Neobalanocarpus heimii (King) Ashton is a highly sought after tree species endemic to the Malay Peninsula. Due to its scarcity and high value, the tree is classified under the IUCN Red List categories of Vulnerable. In this study, we assembled the complete chloroplast (cp) genome of N. heimii using data from high-throughput Illumina sequencing. The Chengal cp genome is 151,191 bp in size and includes two inverted repeat regions of 23,721 bp each, which is separated by a large single copy region of 83,801 bp and a small single copy region of 19,948 bp. A total of 130 genes were predicted, including 37 tRNA, 8 rRNA, and 85 protein-coding genes. Phylogenetic analysis placed N. heimii within the order Malvales.

High cryptic species diversity is revealed by genome-wide polymorphisms in a wild relative of banana, Musa itinerans, and implications for its conservation in subtropical China.

BACKGROUND: Species delimitation is a challenging but essential task in conservation biology. Morphologically similar species are sometimes difficult to recognize even after examination by experienced taxonomists. With the advent of molecular approaches in species delimitation, this hidden diversity has received much recent attention. In addition to DNA barcoding approaches, analytical tools based on the multi-species coalescence model (MSC) have been developed for species delimitation. Musa itinerans is widely distributed in subtropical Asia, and at least six varieties have been documented. However, the number of evolutionarily distinct lineages remains unknown. RESULTS: Using genome resequencing data of five populations (making up four varieties), we examined genome-wide variation and found four varieties that were evolutionary significant units. A Bayesian Phylogenetics and Phylogeography (BP&P) analysis using 123 single copy nuclear genes support three speciation events of M. itinerans varieties with robust posterior speciation probabilities; However, a Bayes factor delimitation of species with genomic data (BFD*) analysis using 1201 unlinked single nucleotide polymorphisms gave decisive support for a five-lineage model. When reconciling divergence time estimates with a speciation time scale, a modified three-lineage model was consistent with that of BP&P, in which the speciation time of two varieties (M. itinerans var. itinerans and M. itinerans var. lechangensis) were dated to 26.2 kya and 10.7 kya, respectively. In contrast, other two varieties (M. itinerans var. chinensis and M. itinerans var. guangdongensis) diverged only 3.8 kya in the Anthropocene; this may be a consequence of genetic drift rather than a speciation event. CONCLUSION: Our results showed that the M. itinerans species complex harbours high cryptic species diversity. We recommend that M. itinerans var. itinerans and M. itinerans var. lechangensis be elevated to subspecies status, and the extremely rare latter subspecies be given priority for conservation. We also recommend that the very recently diverged M. itinerans var. chinensis and M. itinerans var. guangdongensis should be merged under the subspecies M. itinerans var. chinensis. Finally, we speculate that species delimitation of recently diverged lineages may be more effective using genome-wide bi-allelic SNP markers with BFD* than by using unlinked loci and BP&P.

Functional role of circular RNAs in cancer development and progression.

Circular RNAs (circRNAs) are a large class of endogenously expressed non-coding RNAs formed by covalently closed loops through back-splicing. High throughput sequencing technologies have identified thousands of circRNAs with high sequence conservation and cell type specific expression in eukaryotes. CircRNAs play multiple important roles in cellular physiology functioning as miRNA sponges, transcriptional regulators, RBP binding molecules, templates for protein translation, and immune regulators. In a clinical context, circRNAs expression is correlated with patient's clinicopathological features in cancers including breast, liver, gastric, colorectal, and lung cancer. Additionally, distinct properties of circRNAs, such as high stability, exonuclease resistance, and existence in body fluids, show promising role for circRNAs as molecular biomarkers for tumor diagnosis, non-invasive monitoring, prognosis, and therapeutic intervention. Therefore, it is critical to further understand the molecular mechanism underlying circRNAs interaction in tumors and the recent progress of this RNA species in cancer development. In this review, we provide a detailed description of biological functions, molecular role of circRNAs in different cancers, and its potential role as biomarkers in a clinical context.

Comparative transcriptome analysis of the invasive weed Mikania micrantha with its native congeners provides insights into genetic basis underlying successful invasion.

BACKGROUND: Mikania micrantha H.B.K. (Asteraceae) is one of the world's most invasive weeds which has been rapidly expanding in tropical Asia, including China, while its close relative M. cordata, the only Mikania species native to China, shows no harm to the local ecosystems. These two species are very similar in morphology but differ remarkably in several ecological and physiological traits, representing an ideal system for comparative analysis to investigate the genetic basis underlying invasion success. In this study, we performed RNA-sequencing on the invader M. micrantha and its native congener M. cordata in China, to unravel the genetic basis underlying the strong invasiveness of M. micrantha. For a more robust comparison, another non-invasive congener M. cordifolia was also sequenced and compared. RESULTS: A total of 52,179, 55,835, and 52,983 unigenes were obtained for M. micrantha, M. cordata, and M. cordifolia, respectively. Phylogenetic analyses and divergence time dating revealed a relatively recent split between M. micrantha and M. cordata, i.e., approximately 4.81 million years ago (MYA), after their divergence with M. cordifolia (8.70 MYA). Gene ontology classifications, pathway assignments and differential expression analysis revealed higher representation or significant up-regulation of genes associated with photosynthesis, energy metabolism, protein modification and stress response in M. micrantha than in M. cordata or M. cordifolia. Analysis of accelerated evolution and positive selection also suggested the importance of these related genes and processes to the adaptability and invasiveness of M. micrantha. Particularly, most (77 out of 112, i.e. 68.75%) positively selected genes found in M. micrantha could be classified into four groups, i.e., energy acquisition and utilization (10 genes), growth and reproduction (13 genes), protection and repair (34 genes), and signal transduction and expression regulation (20 genes), which may have contributed to the high adaptability of M. micrantha to various new environments and the capability to occupy a wider niche, reflected in its high invasiveness. CONCLUSIONS: We characterized the transcriptomes of the invasive species M. micrantha and its non-invasive congeners, M. cordata and M. cordifolia. A comparison of their transcriptomes provided insights into the genetic basis of the high invasiveness of M. micrantha.